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Abstract 



We give a detailed treatment of the back-reaction effects on the Hawking spec- 
trum in the late-time expansion within the semiclassical approach to the Hawking 
r — [ , radiation. We find that the boundary value problem defining the action of the 

^-p . modes which are regular at the horizon admits in general the presence of caustics. 

We show that for radii less that a certain critical value Tc no caustic occurs for all 
values of the wave number and time and we give a rigorous lower bound on such 
/\f • a critical value. We solve the exact system of non linear equations defining the 

C^ ' motion, by an iterative procedure rigorously convergent at late times. The first two 

terms of such an expansion give the 0{oj/M) correction to the Hawking spectrum. 



1 Introduction 

The semiclassical treatment of the Hawking radiation was introduced by Kraus and 
Wilczek in [H [2] after which several developments followed. The main interest of the 
treatment is to provide a method to compute the back-reaction effect of the radiation 
on the black hole, or in different words a method which keeps into account the conserva- 
tion of energy, an effect which is completely ignored in the the external field treatment 
of the phenomenon [3], HI |5]. The main idea is to replace the free field modes of the 
radiation by the semiclassical wave function of a shell of matter or radiation which con- 
sistently propagates in the gravitational field generated by the back hole and by the shell 
itself. The shell dynamics was studied in detail in many papers (see [SI [D El El El El HD] 
where also a more complete list of references is found) . In the original semiclassical treat- 
ment P, El [11] the spectrum of the Hawking radiation is extracted through the standard 
Fourier analysis of the regular modes. Later such a treatment was related to the tun- 
neling picture; such an approach gave also rise to several proposals and to controversy 
[la [la [H [la [la [la [111 [la EQl Ell EI] and [22] for a vast list of references. 
We think that the mode analysis is still the clearest and safest way to extract the results 
in the semiclassical approach. 

The present paper is devoted to a detailed analysis of the construction of the semiclassical 
modes and their time Fourier transform. The action related to the modes which are regular 
on the horizon is defined through mixed boundary condition, i.e. a condition on the value 
of the conjugate momentum at t = and a condition at time t on the coordinate r. While 
it is easy to prove that the variational problem in which coordinates are given both at 
time and t does not present caustics, i.e. at most one motion satisfies the variational 
problem, we prove that for the above mentioned mixed boundary condition problem in 
general caustics arise i.e. in general more that one trajectory in phase space satisfies the 
mixed boundary conditions. 

Qualitatively this phenomenon is due to the fact that the time to reach the final value ri 
of the radius is an increasing function the mass of the black hole given an initial value 
tq, but such initial value tq through the condition on the initial momentum is also an 
increasing function of the mass of the black hole thus giving rise to two counteracting 
effects. 

On the other hand we prove that if the end point ri, at which the time Fourier analysis 
will be performed, is less that a critical value, caustics to not occur and we give also a 
rigorous lower bound on that critical value. 
In the original paper [IJ it was argued that the semiclassical approximation is expected 



to be valid, for not too large values of ri. If we stay below the critical value r^ we are 
in the favorable situation of absence of caustics where the semiclassical wave function is 
well defined. It is well know on the other hand that the time Fourier analysis gives result 
independent of ri [1] . 

We come now to the computation of the action as a function of t. Such a problem 
corresponds to the solution of a system of two highly non linear equations where the 
two unknown are the value H of the hamiltonian, which even if a constant of motion 
depends on the time t of the boundary condition, and the shell position at time t = 0, 
tq which also depends on the mixed boundary conditions. In P and [9] a discussion 
was given of a simplified system of equations obtained by keeping only the most singular 
terms in the exact equations; even if such a truncation does not remove the high non 
linearity of the system it simplifies it a lot. Instead here we treat the full exact system 
of equations. By introducing an implicit time variable we show that for large times such 
a system of equations is equivalent to an other non linear equation which can be solved 
by a convergent iterative procedure. We show that the first two terms of the convergent 
iterative procedure are sufficient to provide the leading spectrum of the radiation and 
its back-reaction correction terms of order u/M. The treatment is detailed enough to 
show directly how at late times, higher and higher momenta of the regular modes on the 
horizon contribute. From the general viewpoint confining ourselves to large times is not 
a limitation. The reason, as well known, is that the Hawking radiation is a late-time 
phenomenon. This feature was particularly enlightened in the treatment of [5] where 
the Hawking radiation is extracted by the Fourier analysis of the modes on a bounded 
space-time region fixed in space but translated at asymptotically large times. In fact the 
persistence in time is one of the key features of the black hole radiation. 
The paper is organized as follows: sections |2] to S] are devoted to the gauge choice, to 
the description of the reduced action for the shell dynamics and the ensuing equations of 
motion. In section E] we prove the existence of caustics and give a rigorous bound on the 
value of the critical radius Tc below which no caustic develops. Section E] is devoted to the 
non linear system of equations related to the regular modes and to its solution through 
a convergent iterative procedure. Section [7] deals with the saddle point calculation of 
the Bogoliubov coefficients and discusses the region of validity for such a procedure. In 
section IH] we give some concluding remarks. In the Appendix we collect the most important 
formulae relative to the shell dynamics. As it is usual in this field, we work with c = 1 
and G = 1; it means that time, momenta, mass and energy are all measured in units of 
length. 



2 Choice of gauge and the conjugate momentum 

In the general expression of the metric 

ds^ = -N^dt^ + L\dr + N'dtf + R^dn'^ (1) 

all quantities N, L, N^, R are supposed functions only of the radial variable r and t thus 
realizing spherical symmetry. 

The semiclassical approach is best developed in the Painleve-Gullstrand metric charac- 
terized by setting L = 1 in ([1]). Such a metric has the advantage of being non singular at 
the horizon. After fixing L = 1 one has still a gauge choice on R. In presence of a shell of 
matter one cannot choose R = r. One has several choices; for a discussion see [Bl El fTU] . 
In the present paper we will use the "outer gauge" which is defined by i? = r for r > f 
where f denotes the shell position. At r = f, R is continuous as all the other functions 
appearing in ([1]), but its derivative is discontinuous. For the reader's convenience we 
report in the Appendix the main results on the shell dynamics which are necessary in the 
following developments. 

The first step is to go over from the standard Hilbert-Einstein action added to the action 
of the matter shell, to the action expressed in Hamiltonian form. 

As usual in gravity it is better to work on a bounded region of space-time. The radial 
coordinate will range from Tj to r^, while time ranges from tj to tj. 

After solving the constraints one can rigorously express the action in reduced form i.e. 
a form in which only the coordinate f of the shell and a conjugate momentum appears, 
in addition to the boundary terms. As always these are essential in gravity, where the 
boundary terms play the role of the Hamiltonian. The reduced action in the outer gauge 
is given by [H [HI [9] (see also the Appendix) 

S= I ipcf- M{t) I -^dr - HN{r,) + MN{r,) | dt (2) 

where F is the generating function 



F = R^I(Rr-l + ^ (3) 



As a consequence of the constraints the quantity Ai which appears in ([2]) is constant in 
r except at the position of the shell where it is subject to a discontinuity. M and H are 
the the value of the quantity A^ below and above the shell position and thus at r^ and rg. 



One can consider either M or H as a. given datum of the problem. In the outer gauge for 
which the action has the form ([2]) it is simpler to consider M as a datum of the problem 
which to be consistent with the gravitational equations has to be constant in time [S]; so 
the term proportional to M disappears and we reach 

S= r [p^f^- HN{re) + MN{ri))dt. (4) 

Jti 

As in the variation, the components of the metric have to be kept constant at the bound- 
aries, action (j4]) with the normalization N{re) = 1 is equivalent to 

S= / {pcf-H)dt. (5) 

Even though the kinetic momentum of the shell p (see Appendix) is a gauge dependent 
quantity, the conjugate momentum pc appearing in the reduced action is a gauge invariant 
quantity [9]. It can be computed both for a massive or massless shell [H [3 [H [9] ; as we 
shall in this paper be interested in the massless case we report below its expression only 
for the massless case 



Pc = V2Mr - V2Hf - f log -^ ^=. (6) 

One has to keep in mind that pc is not the kinetic momentum of the shell but the conjugate 
momentum with respect to f of the whole system. 

3 The action for the modes regular at the horizon 

As in the following of the paper only f would appear, we will for notational simplicity 
denote f (the shell position) simply by r without any possibility of confusion. 

At the semiclassical level the modes which are invariant under the Killing vector ^ are 
simply given by 

with /p = Gh the square of the Planck length and 

Pcdr — Ht + const. (8) 

As is well know such modes have the feature of being singular at the horizon; this is 
immediately seen from the expression of pc eq. ([6]) which diverges at r = 2H. The vacuum 



given by a^jlO) = 0, being a^ the destruction operator relative to the described modes 
gives rise to a singular description at the horizon, while a free falling observer should not 
experiment any singularity jS]- 

Instead the true vacuum should be described in term of modes which are regular at the 
horizon [HII]. Thus the main object which intervenes in the semiclassical treatment is the 
semiclassical expression of the modes regular at the horizon i.e. the action of the system 
which describes an outgoing shell of matter and has the following boundary conditions 
[1]: i) at time the conjugate momentum is a given value k; ii) at time t the shell position 
r is a given value ri. The expression for such an action was already given by Kraus and 
Wilczek in [T]. With the two conditions Pc{0) = k and r(t) = ri the action is 

S{n,t,k) = kro{ri,t,k)+ f {pcr-H{r{t'),p,{t')))dt' = 



= kro{ri,t,k)+ / p^rdt' - H[ri,t,k]t. (9) 

Jo 

The last equality is due to the fact that H along the motion is a constant despite H 

depends on the boundary conditions as explicitely written, tq denotes the value of r at 

time 0; also such a quantity depends on the imposed boundary conditions. Taking into 

account that r and pc depend both on the final time t and on the running time t', and 

denoting with a dot the derivative with respect to t' one has 



Similarly 



dS , 9ro f\ dr , dr . ^ , .. 

7^ = ^Tr+ ip-Tr + p-Tr)^^ = p- lo 

ori or I Jq ^ or I Ori ' 

dS , 9ro , . ,,, , f\ dr . dr. ,, ^^ , , 

The action ^ has to be computed on the solution of the equation of motion,satisfying 
the described boundary conditions. 



4 The equations of motion 

In the outer gauge [H IS], that we adopt here, the equation of motion for r has the form 



dr 2H , , 

— = 1-W 12 

dt y r ^ ^ 

while dpc/dt can be obtained substituting r(t) in eq.([6]). Eq. flT2|) can be integrated in the 
form 

t = 4iJ'log ^ ~ + ri - ro + 2^/2Hri - 2^2HrQ. (13) 

^- V2H 



The boundary condition at t = gives 



^ - V2H _ [" dH' 



< A; = V2Mro - ^2HrQ - r^ log ^ ^ = / ^= (14) 

V '''0 

where 2M < 2H < rg < ri; eq. fll4p together with eq. fll3p should determine completely the 
motion. However as we mentioned in the introduction we shall find that for sufficiently 
large ri caustics arise, i.e. there exist values of t, ri and k for which the boundary 
conditions are satisfied by more than one motion. 

At t = we have S{ri,0,k) = kri which is regular at the horizon and in virtue of the 
equations of motion S remains regular in the time development. 

5 The occurrence of caustics 

It is very easy to show that the standard variational problem in which r is fixed to tq at 
time and to ri at time t presents no caustics. In fact from 

r^ dr 

Jro 1 _ . m. 



we see that t is an increasing function of H. Thus there is at most one H which satisfies 
the boundary conditions. But the value of r at t and H determine completely the motion 
for an outgoing shell. On the other hand the problem iQ which has mixed boundary 
conditions is more complicated. First we note that being, from eq. (1141) 



f dk\ 1 
— = =>0 16 



we have that at fixed k, if is a single valued function of rg and viceversa tq a single valued 
function of H. Moreover we have 

(dk\ f" dH' f2W 1 /-v^ y'dy ^^ 

{drojH Jm 2ro(l-W^)2V '"o 2 J ^ {I - yf 



with 



Combined with eq. fjTGjl it gives 



2H 

To ' 



fdH\ 1 ^, /"^ y'^dy , , 



To investigate the occurrence of caustics we shall compute the derivative of t with respect 
to To under the constraint of constant k. First we note that from eq.f lT^ 



fdt\ n dr' /2 1 1 /'V^ dz , , 



i^\ -^-^- Pl) 



\dr^ 



0/ H 



dt\ /" dt\ f dt\ fdH\ 1 , ^ ^ , 

+ h^ h:- =-. ^[1-^1^2] (22) 



droJk \droJH \dH J ^^\dro J ^^ 1- 

h-ii-v7.)f_j^,, i^^n^v^r-^. (23) 



and 



Thus 



with 



It is easily seen that 

< Ji < 1, < /2. (24) 

The value of eq. fl22|) for ro = ri, due to the vanishing of I2, is the finite negative value 

(P-) in) = ^-j= < (25) 

with Hi solution of eq.f ll4p in which tq has been substituted with ri. On the other hand 
given a value of k and of tq (which through eq. flT^ gives a value of H with 2M < 2H < tq), 
there will always be ri large enough as to make the product Ji I2 larger that 1; this 
because I2 diverges when the lower integration limit goes to zero. Thus at that point 
eq. fl22j) becomes positive while (125 p still has to hold. Summarizing we found that for a 
given k, for large enough ri the derivative (^2^, when tq moves from ri to 2M changes 
sign, thus vanishing at at least one intermediate point. This implies the occurrence of 
caustics [23] . In fact the vanishing of the derivative ( 122|) at the value Tq implies that there 
will be points Tq and Tq on the right and on the left of Tq which give rise to the same 
value of t. Thus we shall have pairs of distinct motions with the same k which reach ri 
at the same time t. (One can also give numerical examples of such pairs of motions). In 
constructing caustics we took ri large enough. We will show now that for ri < re where 
Tc is a critical value, no caustic arises, for any k. 

Below we give a simple procedure to give a rigorous lower bound on Vc- It is very simple 
to show that for ^/a = 1 both Ji and I2 are equal to 1. Setting 

/i = 1 + Ai, /2 = 1 + A2 (26) 



we will prove that 

Ai + A2 < (27) 

for ri less than a value rj, independently of the value of k. Then being Ji > and I2 > 
we have Ii I2 < 1 and thus {-^)k always negative. 

Thus for ri < r^ there will be no caustic i.e. r^ will constitute a lower bound on Vc- With 
regard to the proof of (ITTI) explicit computation of the integrals gives 



A1 + A2 r- 1 /2M ^, / 2M 

- 2 log 1 



1- ^/a ^_ f2M \ ro \ \ ro 



ro 



21og./^ i= — + ^i — + 21og( l-./^v^) < 



a V '"I l-./m^ ./^\/^ V V ^1 

A/ ri V V ''1 

-l-21og./^ ^= — + ^i — + 21og(l-^/^v^) < 

V ''1 V ''1 



_l_2,ogJ2^ 1^ ^ ^, 2^0, I- jm (28) 

V ri j^ _ /2M /2M_ \ \ ri J 

where in writing the two inequalities we used repeatedly 2M < 2H < tq. The last term in 
eq. fl25]) is a decreasing function of 2M/ri and it is less that zero for 2M/ri = 2/10. Thus 
we do not have caustics for ri < lOM and as a consequence we have rigorously Vc > lOM 
for the critical value Vc- A numerical search of eq.( l22|) gives the wider bound Vc > 24 M. 
In [H [9] the approximate system of equation obtained by retaining in eqs. fll3|14p only 
the singular terms i.e. only the logarithms was considered. Also for this approximate 
system of equations, caustics occur for ri sufficiently large. The occurrence of caustics 
for ri > Tc hints at a failure of the semiclassical approximation when we move too far 
from the horizon as the modes would show a discontinuity in ri at the point where more 
than one trajectory in phase space start contributing. On the other hand we will show in 
sectional that for any given pair (ri, k) even for ri > r^ for t sufficiently large no caustic 
occurs. 

In [1] it was proposed to perform the time Fourier analysis at a point ri not too far from 
the horizon, the reason being that there one should expect the semiclassical approximation 
to be reliable. We showed above that for ri < r^ there are no ambiguities in the definition 
of the action and in addition it is well known the time Fourier transform gives results 
independent of ri, thus we shall work with ri < Vc- 



6 The late-time expansion 

In this section we shall give the solution of the equations for H{t) and ro(t) in the form 
of a series convergent for large times. Large times are also the relevant times for the 
treatment of the Hawking radiation. We recall that 2M < 2H < tq < ri. Then from 
eq. (IT^ for ri fixed, t — )■ +00 implies 

^-V2H ^0. (29) 

Looking now at eq. fll4p we must have in the same limit 



7^-V2M^0 (30) 



and as 2M < 2H < ro we have also H — )■ M. We introduce now the implicit time variable 
T = exp(— ^) which due to the bounds on H, for t — )■ +00 tends to 0. Eq. (TT3|) becomes 



^-V2H ^\ ^H \l2H y2Hj ^ ' 

It will be useful for the following developments to use the notation 



h = V2H, m = V2M, vq = y^, A = y/¥^ - m > (32) 

and set 

h — m = Tch, Vq — m = Tcr (33) 

with ch and cr functions of T do be determined. Eq. (l3ip becomes Fi = 1 with 



J. cr-ch r {A-CRT){A + Am+{2cH + CR)T). 

Fi = -j T^r- exp — ^q-^ 34 

A — Tch 2[m + chT)'' -' 



and eq. (TT^ becomes F2 = k with 



n2,_„cr-ch 



F2 = ~CHT{m + crT) - (m + TcrY log — -. (35) 

We want to express /i as a function of the implicit variable T. 

First we note that for T = the system of the two equations Fi = 1, F2 = k has the 

unique solution 

4 = (e^-l)|, 4 = e^| (36) 



where 



...,(_^(4±M) ,3r, 



and that Fi and F2 in a polydisk around T = 0, ch = c'jj, cr = c^j^ are analytic functions 
of T, ch,cr. We have 



dF, 



dc 



R 



-m^- 



1 d 



0.c°.c° ^R ^H ^R 



f + (38) 



and thus according to the imphcit function theorem [2l], in a neighborhood of c^j, c/j will 
be an analytic function of T and cu- Substituting in Fx we obtain the equation 



Fi(T,ch,cr(T,c^)). (39) 



At T = 0, c// = (Pit we have 



dFi 



dc 



H 



dFi dFidcR E, c%. E 1 , 



-H 



and thus c// will be an analytic function /(T) of T in a neighborhood of T = 0. Recalling 
now the definition of ch we have 

h = m + Tf{T). (41) 

Summing up we found that in a neighborhood of T = from the equations Fi = 1, F2 = k 
eq. fHT]) follows. Eq.l HTj) due to the definition of T = exp(— t/(2/i^)) is still an implicit 
equation. We show now that eq. pTj) can be solved by a convergent iterative procedure. 
Being /(T) analytic with /(O) = c^^ > 0, g(t) = Tf{T) will be in a neighborhood of 
T = 0,T > a non negative function of T of Lipschitz type, i.e. < \g(T2) — g(Ti)\ < 
c\T2 — Ti\ for Ti,T2 belonging to such a neighborhood. Moreover g[0) = 0. For t such 
that 

ce"^ < A (42) 

the r.h.s. of eq. 04ip maps the domain ■m<h<A + m = y/r\ into itself. We start 
the iterative process with ho = m. We must give a bound on \hn+i — hn\. We have for 
hn - hn-i > 

|/in+i - hn\ < ce'^il - e~^ ^^1'^ ) < ce~^ t ?~,""' < -— e"^|/i„ - /i„_i|. 

(43) 
and for /i„ — hn-i < we reach the same result. Thus for t satisfying eq. (H2l) and 



C t __L. 

e 2n < 1 (44) 



mm? 



we have a contraction mapping and according to Banach fixed point theorem [25] eq. fj^T]) 
has one and only one solution given by the convergent sequence /i„. 
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2M 




Figure 1: Time development of 2H(t) and ro(t) 

We give in Fig.l a qualitative graph of the behavior in time of 2H{t) and ro(t). 

We work out now explicitely the first two terms of such an iteration procedure; they will 

be sufficient to give the 0{u/M) corrections to the Hawking distribution. With h^ = m 

and r = e^'/^^*^) = e-^/^^^"^) ^e have 



t 



hi = m + c^jjT, h2 = m + c^jjT + ^c^jjTf + 0{t 



771-^ 



and thus for H(t) 



t 



Hit) = M + V2M4,t + i^icVf + 0{t^) 



(45) 



(46) 



Due to eq. flTT|^ the time dependence of the mode which is regular at the horizon, for fixed 
Ti is 

- I H{t')dt' = const -Mt + 4:MV2Mti + trl + 0{tI) (47) 

with Ti = c^^T i.e. for the semiclassical mode we have 



oiS/l] 



J[q{ri)-Mt+iMV2MTi+tT^+0(T^)]/l 



(48) 



where /p = Gh is the square of the Planck length. Thus S at large times behaves as 
—Mt independently of k. On the other hand the Fourier time analysis of e^^l^p contains 
frequencies which are above and below the value M and this is the well known fact that 
the mode of the system which is regular at the horizon does not represent an eigenvalue 
of the energy as measured by a stationary observer at space infinity. The deviations 



11 



from the value M represent the positive and negative frequency content of the radiation 
mode. One has to keep in mind that the action which appears in fllS]) refers to the whole 
system, which includes both the shell and the core. If we want to analyze the modes of 
the radiation we have to subtract from the exponent the background term —Mt. 

7 The saddle point approximation 

As well known and discussed in [H IH] the Bogoliubov coefficients a^^ and f3^k are given 
by 

a^, = c{t,) j dt e^(^+^^*+-*)/'^ /3^, = c{n) j dt e^(5+Mt-.t)/z|, _ (49) 

As discussed in section [5] we will work with ri < Tc. The above integrals will be computed 
using the saddle point method where /p plays the role of asymptotic parameter [26] . 
From what we derived in the previous section, the exponent appearing in the integrands, 
multiplied by — z/p apart from g(ri) which is constant in time and common to both 
coefficients, are respectively 

2mVi + trl + m^{s + 1)tI ± uot with ri = c^^jT (50) 

where we used the notation of eg. (15^ . For the a^k case (i.e. upper sign) the saddle point 
is given by the value of time t which satisfies 

= -Hit) +M + uj = -rriTi -r^ - sr^ + u (51) 

which being u > has solution for real t and thus at a real value of the exponent in 
eq. fHOj) . On the contrary for the f3uik case (lower sign), the saddle point equation 

= -H{t) + M ~uj = -mri -t^ - st^ - u (52) 



771 



has solution for complex t. At such a value of time the exponent ( 150|) (lower sign) equals 

B = -2m^uj - t{T^ + uj)-{s- l)m^Tl (53) 

The solution of eq. (l52p to second order in u, which is the order we are interested in, is 
given by 

^1 = -- i-z:iM--rr-) + -i • (54) 



From eq. flS^ we see that to find the imaginary part of such exponent to order uP' we 
simply need the imaginary part of t to first order in u. Using (I54p we have 

Im t = -27Tm\l - ^). (55) 



m^' 
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Substituting into eg. (153!) we find 



ImB = 2TTm'^uj{l - ^) = 4:7tMuj{1 - -^) (56) 



which according to (H^ has to be divided by /p. Thus we have 



Wujk? 



(57) 



which is independent of k. We see from eq. fH6|) that for t — )■ +oo, H{t) tends to M 
and thus the time Fourier transform of the exponential of the action f l48|) which refers to 
the whole system has a singularity at the frequency M. Recalling that H (outer mass) 
represents the energy of the whole system, we identify the parameter M with the mass of 
the black-hole before the decay. 
Using the property of the Bogoliubov coefficients 

^{aujkoCo'k - l^u^kl^l'k) = ^^,u>' (58) 

k 

one reaches for the flux of the Hawking radiation [11] 



^M^^ = |^ 8.M../.. • (59) 

e p — 1 

This completes the explicit derivation of the u"^ correction to the Hawking formula from 
the time Fourier transform of the semiclassical modes. 



An alternative way to derive (1571) was given by Keski-Vakkuri and Kraus [llj where it is 
proven that for the /3^fc coefficient the imaginary part of the action at the saddle point 
(1521) is given by 

rri /•2M -i 

Im / pjr = Im / pjr = 7i-{{2Mf - {2Hf) = 4:7iMuj{l - — -) (60) 

Jro J2H 2 2M 

which is equivalent to eq. (l56p . The importance of equation (16 Op is to show directly how 
the "tunneling" is due only to the imaginary part of the "space part" of the action. 
With regard to the validity of the expansion we see from the saddle point value (I51f52p 



k 

A(e2M — 1) t u , , 

V2ME 2M ^ ' 

that the series if effectively an expansion in u/M and thus expected to hold for oj/M << 1. 



From eq. fl6T|) we see that for a given cj, large values of the wave number k contribute at 
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times t which grow hke 2k. The typical u for the radiation emitted by a black hole of 
mass M is according to eq. flFTI) (Wien's law) 

u « — ^ (62) 

and thus the approximation expected to be reliable at the typical frequency (162|) or below 
for Ip/SirM'^ « 1 i.e. for black holes of mass of a few Planck masses or of higher mass. 



8 Conclusions 

In this paper we gave a detailed treatment of the late-time expansion which occurs in the 
semiclassical approach to the Hawking radiation. We find that the variational problem 
defining the action related to the modes which are regular at the horizon allows in general 
more than one solution, due to the presence of caustics. We prove however that for radii 
below a critical value Vc the variational problem has only one solution and we give a rig- 
orous lower bound on r^. Thus for ri less that r^ where the semiclassical approximation is 
expected to be accurate there are non ambiguities in computing the action and the time 
Fourier transform can be applied to extract the Bogoliubov coefficients. The Hamiltonian 
depends on the boundary condition through a system of two highly non linear equations. 
We show that for sufficiently late times such a system of equation is rigorously equivalent 
to an other non linear equation which can be solved through a convergent iterative proce- 
dure. We work out explicitely the first two steps of such iteration which are sufficient to 
compute the u/M correction to the Hawking spectrum. The treatment shows directly the 
relation between late times and high wave numbers of the modes regular at the horizon. 
The first two terms in the iterative process are sufficient to give accurate results for the 
back-reaction effects for frequencies at or below the typical frequency of the spectrum and 
black holes of a few Planck masses or higher mass. 
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Appendix 

We summarize here the essential formulas of the shell dynamics. For more details see 
[H [Zl El E] . One starts from the usual Hilbert-Einstein action to which the shell action is 
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added 

S = j^-^ / Ry/^ d^x + S shell- (63) 

We shall in the following use c = G = 1 which simply means that masses acquire the 
dimension of length i.e. they are measured by the related Schwarzschild radius divided by 
2. As usual in gravity it is better to work on a bounded region of space-time. Employing 
the general spherically symmetric metric ([T]) the action can be rewritten in Hamiltonian 
form as P, [I] 

/•*/ r<^ . . /•*/ NRR' 

S = / dt dr{TTLL + ttrR - NHt - N'Hr) + dt {-N'ttlL + ——\ 



dt pf (64) 

where f denotes the radial coordinate of the shell. The constraints are given by 

nr = nnR'-n'LL-p6ir-f), (65) 

?/2 r^2 



RR" R'^ Lnl RR'L' ttlttr L 



The Painleve-GuUstrand gauge is defined by L = 1. There is still one gauge freedom in 

the choice of R{r). In virtue of the constraints R'{r) has to be discontinuous at r = f. 

Here we will adopt the "outer gauge" [S] defined by R{r) = r for r > f i.e. in the massless 

case 

R(r) = r + -g(r — r) (67) 

r 

with g smooth function of support [— /,0], g{Q) = and g'{0—) = 1. Other gauges could 
well be used [H |9] . The constraints can be solved and the action in the outer gauge takes 

the form 

ftf / rf{t) Qp 



S= I ( p, r - M{t) I -^dr - HN{r,) + MN{r,) | dt (68) 

where F is the generating function 

F = RW + RR\C - B) (69) 

with 



iy = ^i?--i + ^, £ = iog(i?'-iy), ^ = V^ + iog(i-V^)- (70) 

The general expression of the conjugate momentum p^ is [S] 

Pc = R{AC - AB) (71) 
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where A represents the discontinuity of the related quantities across the shell position f. 
Contrary to p, p^ is a gauge invariant quantity within the Painleve class of gauges [9]. 
Its expression for the case of a niassless shell is given by eq.([n]). Normalizing the lapse 
function N, which is constant for r > f, as N{re) = 1 we have from the expression ([6]) of 
Pc and action ([5]) the equation of motion 



dH ^ 


l2H , 




\ = T 


dpc 


V f 



(72) 
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